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We study the dynamical generation of entanglement for a quantum system consisting in a pair of 
interacting spins, si and S2, in a constant magnetic field. Three different situations are considered: 
a) si = S2 = 1/2, b) si — > oo,S2 — 1/2 and c) si = S2 ^ oo, corresponding respectively to a fully 
quantum system (two interacting qubits), a quantum degree of freedom coupled to a semiclassical 
one (a qubit in contact with an environment) and a fully semiclassical system, which displays 
chaotic behavior. The time evolution of entanglement for initially separable states is measured 
using the von-Neumann entropy and the linear entropy, and a comparison is made between the two. 
Compactness of the phase space and chaos in the classical limit are seen to have a strong influence 
on the recoherences of the system. 

PACS numbers: 03.65.Ud, 05.45.Mt 



Recent technological advances have made it possible to create and manipulate individual quantum states in the 
laboratory, thus allowing direct observations of entanglement and decoherence [????????], concepts that are 
central to quantum computation and quantum information [? ? ] . Two subsystems 1 and 2 are said to be entangled 
if it is not possible to write the total system density matrix in the form 

P^J^P'Pi^P^^ (1) 

i 

and decoherence is the suppression of off-diagonal elements in a subsystem reduced density matrix [? ] . Quantifying 
entanglement, however, is a difhcult task. Many different measures have been proposed, such as entropy of formation, 
entropy of distillation, and many different distances [?????]. All of them are very hard to calculate in practice, 
although Wootters [? ] has provided a closed form for the entropy of formation in the case of a system consisting 
of only two qubits. If the initial state is separable, the entropy of formation for a bipartite system reduces to the 
von-Neumann entropy (we reserve the letter S for spins) 

= -Tr [p2 logrf p2] , (2) 

where p2 is the reduced density matrix of the subsystem 2 and d is the dimension of its Hilbert space. Another 
easy-to-handle measure is the linear entropy (or idempotency defect) 

6=^il-TT[ip2f]). (3) 

When one considers only globally pure states (as we do here), both the linear and the von-Neumann entropy are 
for separable states and 1 for maximally entangled states. If the initial state is separable, one is in fact considering 
a dynamical, deterministic, generation of entanglement. In this situation how fast entropy grows indicates how fast 
one subsystem suffers decoherence due to entanglement with the other. 

Another interesting question is about the dependence of quantum entanglement upon the classical dynamics of the 
system, more specifically what is the role of chaos in the decoherence process. This has been the subject of many 
recent investigations. In a seminal work Furuya et al [? ] have studied the Jaynes-Cummings model without the 
rotating-wave approximation, and found the entanglement rate to be greater for chaotic initial conditions. In later 
works [? ? ] they have shown that a regular initial condition can sometimes lead to faster entanglement than a 
chaotic one and that recoherences were related to the shape of the spin orbits. In the past few years much work has 
been done in this area, considering continuous variables [? ? ], tops [????], and spin chains [? ? ? ]. 

In the present work we consider two spins. Si and ^2, in a constant magnetic field B = B^z, interacting according 
to the Hamiltonian 

H = eiBnSl + e2BoS'2+aSlS^, (4) 

which is somewhat reminiscent of the usual XXZ Hamiltonian widely used in statistical mechanics. Spin Hamiltonians 
have been considered in connection to quantum computation proposals using NMR [? ? ] (see also [? ]) and cold 
atoms in optical lattices [???]. Entanglement in spin chains at finite temperature was considered for example in [? 
? ]. For simplicity we shall assume ei-Bo = £2^0 — 1- This very simple case already shows a rich behavior depending 
on the parameter a and the spin magnitudes si and 82- 
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The article is organized as follows: in section I we consider si = S2 = 1/2, which corresponds to a pair of qubits. 
In section II we let Si — *■ oo to simulate an environment and study its ability to induce decoherence upon spin 82- 
Section 3 deals with the relation between chaos and decoherence and we conclude in section IV. 



I. TWO QUBITS (si 



S2 



1/2) 



Let us first consider the case si = S2 — 1/2. The Hamiltonian in this case is a 4 x 4 matrix that can easily be 
diagonalized. Its eigenvectors, written in the standard st := {|11), |10), |01), |00)} basis are 
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with respective eigenvalues 



where 



ei = ^, 62 = --, 63 = p, 64 = -p, 



(3 = 



Va2 + 16 



(6) 



(7) 



In order to study dynamical entanglement generation, we use several different initially separable configurations. At 
t = 0, the spins si and S2 are in the states 



|0) + z.|l) 



i = 1,2 



and, after some symmetry considerations, we have decided to investigate the following values of the parameters: 

a)zi — Z2 — 0, b)zi — 1,Z2 = 

c)zi — > 00, Z2 = 0, d)zi = 0, Z2 = 1 

e)zi = 0, Z2 = i, f)zi = Z2 = l 
5)2:1 = 1,22 = 1, h)zi = Z2 = i 



(8) 



(9) 



Note that zi = ±1 and zi = ±i correspond to eigenstates of Sf and Sf, respectively, and the same holds for spin S2- 
A measure of entanglement exists for two spins: the entropy of formation, which measures "the amount of entangle- 
ment necessary to create the entangled state" [? ]. If the system is initially separable, which is the case we consider 
here, the entropy of formation is equal to the von-Ncumann entropy. For each one of the chosen initial conditions 
we followed the time evolution of this quantity and of linear entropy. Both measures have approximately the same 
form, with maxima and minima at the same points, even tough they are not proportional to each other. We have also 
observed that 6 < 6n for all chosen initial conditions, but this is not a general property, since the opposite occurs in 
section III. In Fig. la we see the case a), which is very similar to cases b) — e), all in which at least one of the spins 
is in the |0) state. Conditions /) — h) are also quite similar, all in which \zi\ = [22 1 = 1- Condition /) can be seen in 
Fig. lb. All time scales present in the system can be related to the eigenvalues and are therefore determined by 
a. 



II. A QUBIT AND AN ENVIRONMENT (si ^ oo,S2 = 1/2) 

We now consider S2 = 1/2 and si — 200. Such a huge value of spin is not intended to be realistic, but only to 
simulate an environment for dissipation, when infinitely many degrees of freedom are required (see [? ] and references 
therein). We consider spin si to be, at i = 0, in three different conditions: in a coherent state [? ] 

(10) 
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in a non-localized state 

|V) = (2si + l)-i £ K) (11) 

nil— — Si 

and in a thermal mixture 

Pi = ^ E e—^/'^KX™!!, (12) 

mi=-si 

where iV is a normalization factor (the partition function) 

Let us first consider both spins to be initially in coherent states with zi = Z2 = 0. The linear entropy is shown in 
Fig. 2 (the von-Neumann entropy is indistinguishable from it in this scale) . We see that the entanglement reaches its 
maximum value very fast, and remains saturated for some time, but a strong and localized recoherence takes place in 
a nearly periodic way. As already observed in [? ] , such recoherences are related to the mean value of the z-component 
of the spin, (S"!), or equivalently to how close is the spin to the north pole of the Bloch sphere [? ]. In other words, 
recoherence is an effect of the compactness of the phase space. 

In Fig. 3a we see the results for zi — 0, Z2 = 1. Even though the value of Zi is the same as in the previous case, the 
entanglement is completely different. The linear entropy is now less than the von-Neumann entropy, but they have 
essentially the same information. Entanglement increases with time, but much slower than the previous case. Instead 
of raising immediately to its maximum value, it increases by jumps, alternating fast increments with plateaus. 

We also consider = Z2 = 1, shown in Fig. 3b. The behavior of the entropies is quite similar to the previous case, 
but in a very different scale: entanglement is nearly a hundred times smaller in this case, i.e., the system is practically 
entanglement-free. It is clear that the entanglement rate strongly depends on the initial state of the environment, 
even if this is a coherent state. 

Finally, we have considered the environment to be initially in the non-localized state Hll(l and in the thermal mixture 
(|12|l . In all cases we have used Z2 = 1, and in Fig. 4 we show the behavior of the linear entropy. The qualitative 
behavior is the same for all initial states (in the thermal mixture we have used T = si/10), but the ability to induce 
decoherence is considerably greater for the coherent state. 



III. ENTANGLEMENT AND CHAOS (si = sa ^ oo) 

We now consider both spins to be large (semiclassical), namely we use si = §2 = 15. For any spin Hamiltonian it is 
possible to obtain a semiclassical dynamics using su{2) coherent states (|10|) . in which the quantum dynamics is then 
approximated by Hamilton equations of motion [? ] 

where H — {z\H\z) plays the part of a classical Hamiltonian. In our case this is given by [? ] 

1 



n=-iAi+A2- 31) + -9ig2v/(4si - Ai)(4s2 - A2), (14) 



where 



and the canonical coordinates are defined by 



A=qf+pl (15) 



1,2. (16) 



12' 



In Fig. 5 we see a Poincare section showing the (gi,pi) coordinates at p2 = 0. We use the marked points as initial 
conditions for our coherent states to be evolved. The trajectory marked with a triangle is chaotic, the one marked 
with a square is regular and the circle denotes a periodic orbit. Note that we do not evolve the state semiclassically: 
a semiclassical evolution of the kind H13|l would preserve, by construction, the coherence of the initial state, and thus 
produce no entanglement. 
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In Fig. 6 wc sec the linear entropy as a function of time. The rate of entanglement for short times is very similar 
for all initial conditions, contrary to what is observed for example in the Jaynes-Cummings system [? ]. We also note 
that the maximal value of S is higher for the periodic orbit than for the more generic regular trajectory. Numerical 
analysis indicate that both the short-time entanglement rate and the maximum entanglement value are related not 
to the classical dynamics but to strictly quantum properties of the wave packet, such as energy uncertainty. 

There is however a strong difference between the regular and chaotic initial conditions: the former present strong 
recoherences that are not present in the latter. As already observed in the previous section these recoherences are 
related to the mean value of the 2;-component of the spin, but in this case the relevant spin is si, as we see in Fig. 
7. Note that coherent states have localized Husimi distributions in phase space and thus some more general, less 
localized state, would be less likely to display recoherences. 

Again the entanglement measures are qualitatively very similar, even tough they are not proportional, but in this 
case the von-Neumann entropy is less than the linear entropy at all times and for all initial conditions. We show this 
only for the periodic orbit in Fig. 8. 

IV. SUMMARY 

We have studied dynamical entanglement generation for a pair of interacting spins in a magnetic field. The case of 
two qubits is rather simple, since it allows an analytic treatment up to a point and is useful as a background against 
which we can compare the other cases. We noted for example that for all initial conditions S < 6n- The case of a 
qubit interacting with an environment had a much more rich variety. We have seen that the decohering power of 
the environment strongly depends on its initial state. A comparison was made between different coherent states and 
between a coherent state, a non-localized state and a thermal mixture. 

In the case of two semiclassical spins we have analyzed the influence of chaos upon the entanglement process. We 
saw that chaotic initial conditions presented no recoherences, while regular ones do so. We also noted that in this 
case S > 6n for all initial conditions, contrary to the low dimensional case of two qubits. In all cases the recoherences 
present in time evolution were related to the compactness of spin phase space and localization of coherent states. It 
would bo interesting to investigate other measures of entanglement, specially for the case of initially mixed systems. 
Such work is underway. 
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FIG. 1: (si = S2 = 1/2) Linear (dashed line) and von-Neumann (solid line) entropies for zi = Z2 = (upper panel) and for 
zi = Z2 = 1 (lower panel). We have observed that Sn > S for all times and all initial conditions. 
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FIG. 2: (si = 200, S2 = 1/2) Linear entropy (solid line) and angular momentum (dsished line) as functions of time for 
zi = Z2 = 0. We see that S saturates very fast, but strong recoherences appear, related to (S^)- 



7 




to 
m5 




0,005 



0,000 



FIG. 3: (si = 200, ,S2 = 1/2) Linear (dashed line) and von-Neumann (solid line) entropies for zi = 0, z-z = 1 (upper panel) and 
for zi = Z2 = 1 (lower panel). In both cases we see plateaus followed by jumps, but there is a large difference in magnitude. 
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FIG. 5: Poincare section of the classical Hamiltonian. We use the marked points as initial conditions for our coherent states 
to be evolved. 
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FIG. 6: (si = S2 = 15) Linear entropy 5{t) for three different initial conditions. The regular cases present strong recoherences, 
while these are suppressed in the chaotic case. 
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FIG. 7: (si = S2 = 15) Linear entropy and angular momentum ai = for the periodic initial condition. Recoherences 

are related to the compactness of phase space. 
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FIG. 8: (si = S2 = 15) Linear (dashed line) and von-Neumann (solid line) entropies for the periodic initial condition, 
contrast to section I, now 5 > 5n for all times and all initial conditions. 



